Band structure of semimagnetic Hgi_ y MuyTe quantum wells 
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The band structure of semimagnetic Hgi-yMnyTe/Hgi-^Cd^Te type-III quantum wells has been 
calculated using eight-band k ■ p model in an envelope function approach. Details of the band 
structure calculations are given for the Mn free case (y = 0). A mean field approach is used to 
D ' take the influence of the sp — d exchange interaction on the band structure of QW's with low Mn 

^0 \ concentrations into account. The calculated Landau level fan diagram and the density of states of a 

Hgo.98Mno.02Te/Hgo.3Cdo.7Te QW are in good agreement with recent experimental transport obser- 
vations. The model can be used to interpret the mutual influence of the two-dimensional confinement 
and the sp — d exchange interaction on the transport properties of Hgi-^MnyTe/Hgi-zCdzTe QW's. 



PACS numbers: 71.20.Nr, 71.70.Ej, 71.70.Gm 

'd 

C/2 ■ I. INTRODUCTION 

<D 

Recently, numerous spin related observations have been published which involve optical and transport experiments 
on diluted magnetic, semiconducting heterostructures and quantum wells (QW's)4*2*&^£ The correct interpretation 
of these effects requires a detailed knowledge of the underlying band structure. This is especially important for narrow 
gap semiconductors^ because strong band mixing prevents a simple interpretation of optical and transport results 
1 , ■ by means of a parabolic band model which might be still applicable for most wide gap materials such as GaAs or 
InGaAsAa 

Here, we concentrate on the band structure calculations of HgTe/Hg 1 _ a; Cd a; Te quantum wells. This material has 
some interesting properties: depending on the QW width (dw), the QW has either a normal or inverted band structure 
when dw < 6 nm or dw > 6 nm, respectively. In the latter case the conduction band exhibits Tg symmetry which 
leads to a strong Rashba spin-orbit splitting in QW's with an asymmetrical confinement potential**^ Additionally, 
, the spin splitting of the subbands can be enhanced by introducing magnetic ions (Mn) in the QW structure, e.g., 
■ Hgi_ y Mn y Te /Hgi^^Cd^Te . It should be noted that in II- VI semiconductors, Mn is incorporated into the crystal 
ON I lattice isoelectrically and does not act as a donor or an acceptor. Therefore, Mn ions act primarily as a magnetic 
but not as a Coulomb impurity and mobilities achieved for these QW structures with low Mn concentrations are 
comparable with those for non-magnetic structures. 

This paper is organized as follows: In section II a detailed description of the model used for the band struc- 
ture calculations is presented. In subsection II-A, we consider the model for non-magnetic as well as for mag- 
netic QW's at zero external magnetic field. This model is used to calculate the subband energy dispersion of 
Hgo.98Mno.02Te /Hgo.3Cdo.7Te and HgTe/Hg .3Cdo.7Te QW's. In subsection II-B, the band structure model for 
a HgTe/Hgi-^CdzTe QW in an external magnetic field is described. This model is extended in subsection II-C in 
1 ■ order to take the influence of the sp—d exchange interaction on the band structure of magnetic QW's into account. 
The Landau level fan diagram and the density of states of Hgo.9sMno.02Te /Hgo.3Cdo.7Te QW are compared with 
recent transport experiments, and section III summarizes the results. In the appendices details of the calculations for 
different growth directions (a) and strain as well as piezoelectric effects (b) are discussed. 



' Also at School of Physics and Astronomy, University of Nottingham, Nottingham NG7 2RD, UK; and Physics Department, University 
of Texas at Austin, Austin, Texas 78712-0264, USA. 
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II. BAND STRUCTURE MODEL 

A. 5 = 

The band structure model we use is based on an envelope function approximation introduced by BurtiiS The total 
wave function is expanded in terms of band edge (k = 0) Bloch functions u n : 



#(r)=^F n (r) u n (r), 



(1) 



where F n (r) are the envelope functions. u n is assumed to be the same in the barrier and the well layers. Assuming 
translation invariance in the plane perpendicular to the growth direction [z axis) F n can be represented as follows 



F n = exp[i(k x x + k y y)} f„(z), 



(2) 



where k x and k y are the wave vector components in the plane of the QW. The envelope functions and the energy 
levels near k = are determined within the framework of k ■ p theory by solving a system of coupled differential 
equations ! 11 ^ 2 



E E k <* D nn>fy + E P nn'k a + E n , ( Z )S nn > + V{z)5 nn , ) /„,(*) = E ■ f n (z), 
n f V a, {3 a 



(3) 



where n and n' are the summation indices for the sum over the dimensionality of the chosen basis set, E n i (z) are 
the respective band edge potentials, and V{z) is the self-consistently calculated Hartree-potential. The momentum 
matrix elements P" n t, which describe the coupling between the n and n' bands, are treated exactly, while the D"^,, 
which describe the perturbative coupling of the n and n' bands to the remote bands, are calculated using second-order 
perturbation theoryiiii^ 

In narrow-gap HgTe based structures, the strong coupling between s-like conduction and p-like valence bands causes 
mixing of the electronic states and induces nonparabolicity in the conduction bands. These effects were taken into 
account exactly by Kanei^ in the framework of the k ■ p theory. In order to consider the coupling between the Tg , Tj 
and Tg bands we choose the usual eight-band basis set (see Refs. Eland EH : 

ST 
si 

l/V2(X + iY) t 
l/V6[{X + iY) i -2Z t] 

-l/V&[{X -iY)] +2Z [] (4) 
-l/V2(X-tY) i 
l/V3[{X + iY) i +Z T] 
1/V3[(X - iY) 1 -Z I]. 

The total angular momentum is then given by j — 1/2 or j — 3/2. 

For the chosen basis set, the Hamiltonian H nn > in Eq. [21 for a two-dimensional system with [001] growth direction 



tii(r) = 


|r 6 , 


+ 1/2) 


u 2 {r) = 


|r 6 , 


-1/2) 
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|r 8 , 


+3/2) 


u 4 (r) = 
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-1/2) 
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takes the following form:' 
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T 


Pk z 








T 


3 Pkz 

~7s Pk - 
7s P ^ 



where 





u + v 

-si 





-7E Pk + 

u-v 

c* 
i?t 

V2V 
/Jet 



7e Pk ~ 

\fi Pkz 

R 

C 

U-V 

s+ 

V 2 J + 






1? 



5 



~^ Pk * 
-^ Pk + 

V2V 



u + v 

V2R 



i 



4^5 



+ 



U-A 



-V2R 
-V2V 

c 

U-A J 



1.2 

II 

T 
U 
V 
R 
S± 

s± 

c 



k± — k x i iky, k z 



-id/dz, 



E c {z) + A- ((2F + l)fe 2 + k z (2F + l)k z ) , 



E v (z) - ^— ^7ifc|| + k z ~fik z j 



2mo 

h 2 

2mo 

h 2 

2rriQ 

h 2 



(j2 k f\ - 2k z j 2 k z 

VSfikl - Vs^k 2 ^ , 
V%k± ({73, k z } + [k, k z }) , 



., ,, lo ^ k ± K73,M - k z ]j , 



ft 2 
m 



(5) 



(6) 



[A, B] — AB — BA is the usual commutator and {^4, B} = AB + BA is the anti-commutator for the operators A and B; 
P is the Kane momentum matrix element; E c (z) and E v (z) are the conduction and valence band edges, respectively; 
A is the spin-orbit splitting energy; and 71, 72, 73, n and F describe the coupling to the remote bands and result 
in the fj, and 7 parameters according to /i = 1/2 (73 — 72) and 7=1/2 (73 + 72)- Only the terms with non-spherical 
(cubic) symmetry in the Hamiltonian are proportional to the warping parameter fj,. The case of /i — corresponds 
to the axial approximation. The intrinsic inversion asymmetry is neglected in the Hamiltonian because this effect is 
very small in HgTe based structures^ The band structure parameters for HgTe and CdTe at T = K are listed in 
TableQ] The dependence of the band gap (E g ) of Hgi-^Cd^Te on the temperature and composition x is determined 
from the empirical expression according to Laurenti et alJ& The valence band offset between HgTe and CdTe is taken 
to be equal to 570 meV at T — K, in agreement with recent experiments, 18 and is assumed to vary linearly with 

^.19 



TABLE I: Band structure parameters of HgTe and CdTe at T = 





E g 


A 




F 


71 


72 


73 




HgTe 


-0.303 eV 


1.08 eV 


18.8 eV 





4.1 


0.5 


1.3 


-0.4 


CdTe 


1.606 eV 


0.91 eV 


18.8 eV 


-0.09 


1.47 


-0.28 


0.03 


-1.31 



So far, we have only considered the case of HgTe/Hgi^^Cd^Te QW's with (001) orientation. However, HgTe based 
structures have also been investigated with orientations other than (001), for example, (112) heterostructurespiSi 2 ^ 
The electronic properties of such systems depend strongly on the growth direction. An extension of the model to QW's 
of a given (kkl) orientation can be obtained using the approach of Los et a/i2I: the set of basis functions [Eqs. |@}] 
is changed to a set which adopts the symmetry of the problem, and thus the transformed Hamiltonian once again 
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has the form of Eq. while the matrix elements [Eqs. 10] contain additional terms depending on the structure 
orientation^ The exact formulas for these terms are given in appendix A. Also strain and piezoelectric effects can be 
included, as discussed in appendix B. 

During the last two decades much attention has been paid to the theoretical and experimental understanding of 
diluted magnetic semiconductors (DMS), both in bulk 22 23-24 as well as in low-dimensional structures iS 4 * 2 ^ Extensively 
studied examples of this category are A I 1 I _ J/ Mn 9 B VI alloys, in which the group II component is replaced randomly by the 
transition metal Mn. 26 So far, most research on magnetic two-dimensional structures has been done on wide gap DMS 
materialsiSii In the present work we consider QW's with magnetic ions (Mn) in narrow gap Hgi_ y Mnj,Te /Hgi-^Cd^Te 
QW's. The two-dimensional confinement in the DMS based layer combined with the exchange interaction between 
localized Mn magnetic moments and mobile band electrons make such structures quite interesting candidates for the 
study of their electronic and magnetic properties. 

The band structure of Hgi-yMn^Te /Hgi-^Cd^Te QW's in the absence of a magnetic field can be calculated 
similarly as that of non-magnetic QW's, cf. the description above. The only difference is that the band structure 
parameters for the well now depend on the Mn concentration^ 

In Fig. n the zero-field subband dispersion E(ki\) of n-type Hgo.9sMno.02Te/Hgo.3Cdo.7Te (a, b, d) and 
HgTe/Hgo.3Cdo.7Te (c) (001) QW's are presented. The temperature T and the QW's width dw are set at 4.2 K 
and 12.2 nm, respectively. Fig.Q](a) corresponds to a Hgo.9sMno.02Te/Hgo.3Cdo.7Te QW with symmetrically n-type 
modulation doped barriers, while only the barrier on the substrate side is doped for the case presented in Fig. ^(b). 
In the calculations it is assumed that (i) all donors in the doped layer are ionized, (ii) the charge density is constant 
in this region, and (iii) all electrons are transferred to the QW. The depletion charge in the doped layers is taken to 
be rtDL = — T-2DEG for an asymmetrically and bdl = — 1/2 ro2DEG for a symmetrically doped structure, respectively^ 
("-2DEG denotes the charge density of the 2DEG in the QW and is chosen to be 712DEG = 3.47 x 10 12 cm~ 2 for the 
calculations presented here.) The eigenvalue [Eq. and Poisson equations for the two-dimensional charge carriers 
in the QW are solved self-consistently for both cases. The depletion charge in the doped layers, which is assumed to 
be fixed at the levels indicated above, is included into the boundary conditions to solve the Poisson equation^ Then 
the Hartree potential is determined according to the charge distribution of electrons in the QW, which is given by 
the summation over all conduction band states i and all components n of the envelope functions f n {z): 



where e is the electron charge, /f(E) is the Fermi function. Here, we assume that tt-2DEG < and tidl > 0. In 
analogy, we have for the holes: 



In this case the summation index i runs over all valence band states, and tidl < 0. The self-consistent Hartree 
potential for zero magnetic field is then used to calculate the Landau levels of the 2DEG (see description below). 

From Fig. ^ ( a ) and (b) one observes that both QW's exhibit an inverted band structure. The conduction band 
includes two occupied subbands (labeled as HI and E2); the lower conduction subband (HI) exhibits a heavy hole 
character at fen = 0. The electron-like El subband lies, in this case, below the H2 subband and is now one of the valence 
subbands. The HI and H2 subbands are not split for the symmetric QW [case (a)]; however, for the asymmetric QW, 
a spin-orbit splitting (in the following denoted as a spin splitting as is common in the literature) of the H2 subband 
as well as a pronounced splitting of the HI subband are visible. The splitting of HI is 33.4 meV at kpi and 30.9 meV 
at kpi 1 respectively. The carrier concentrations in the HI— and H1+ subbands for the asymmetric QW, after self- 
consistency has been obtained, are 1.18 x 10 12 and 1.55 x 10 12 cm -2 , respectively. This large spin-orbit splitting, 
usually called Rashba splitting, of the HI state in an asymmetric QW was demonstrated to be an unique feature of 
type-Ill QW's in the inverted band regime^ Experimentally values up to 30 meV have indeed been observed^ 

The subband dispersion E(k\\) of HgTe/Hgo.3Cdo.7Te and Hgo.9sMno.02Te /Hgo.3Cdo.7Te QW's are shown in Fig.Q] 

(c) and (d), where 7Zdl is taken to be equal to — 0.34-n2DEG- The band structures of the non-magnetic (c) and magnetic 

(d) QW's differ notably in the subband separation at fey = 0; the separation between the E2 and HI subbands is about 
24 meV smaller for non-magnetic QW. The Rashba splitting of the HI subband is 15.5 meV for the non-magnetic 
and 13.1 meV for the magnetic QW's for fey (1,0). The corresponding values are 14.2 meV and 12.2 meV for fey(l, 1). 
This relatively small difference arises because the actual band gap (Ajji-m) changes only by 2% upon introduction 
of Mn. 




(7) 




(8) 
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FIG. 1: Band structure of n-type (a) symmetrical and (b,d) asymmetrical Hgo.98Mno.02Te ZHgo.3Cdo.7Te QW's, as well as 
(c) an asymmetrical HgTe/Hgo. 3 Cdo. 7 Te QW, all with d w = 12.2 nm at T = 4.2 K. k F1 = k F (Hl-) and k F2 = k F (Hl+). 
Here the fcy vector is fey (1,0), however the difference between E(ku(l, 0)) and E(k«(l,l)) is less then 6 meV at k F and the 
fc-dependence shows qualitatively the same behavior. The El subbands are shown as dashed lines. 



B. B ^ 

In an external magnetic field perpendicular to the plane of the 2DEG, the electronic bands are split into a series of 
Landau levels. The effects of a magnetic field B = (0, 0, B) can be incorporated by a Peierls substitutioni^SS in the 
Hamiltonian [Eq. 10] as follows: 

k -» k! = -iV + e -A : (9) 
h 

where A is the magnetic vector potential, B = V x A. One possible Landau gauge for B\\z is A = (0, Bx,0). The 
operator k' satisfies the following gauge-invariant relation: 

k' x k' = -i-B. (10) 
h 

From now on, we drop the index ' for simplicity. Subsequently, k x and k y are rewritten such thari^i^ 

a=-j=-k-, a t = -^=-fc+, (11) 

where l c = \J~-^ is the magnetic length; and a and are, respectively, the annihilation and creation operators for 
the harmonic oscillator functions ip n , wherein 

atp n = \/ntp n _ 1 , 

a) (p n = Vn + ltfn+l, (12) 

a^aifn = mp n . 



6 



Here, n = 0,1,2,... are the eigenvalues of the operator a^a. Thus, we can present the Hamiltonian in Eq. (JSJ as a 
function of a, a\ k z = —id/dz, the band structure parameters and their z-dependence. 

Additionally, the Zeeman term, H z , has to be included in the Hamiltonian [Eq. ©]. As shown by Weiler^ this 
leads to the following matrix: 
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(13) 
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In the axial approximation we now assume that the total wave function can be written as& 



/ fi(z) <fm \ 

h{z) <Pn 3 

Pj) h{z) <Pn 5 

h{z) <Pn 7 

V h{z) Vn 8 J 



*Af(r) = CXp 



= exp 
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AN) 
■16 
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AN) 
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[z) 


<Pn 


\ fr 


[z) 





(14) 



where restrictions on the quantum numbers N on the right-hand-side can be derived straightforwardly from Eqs. I|12|l . 
Since n = 0, 1, 2, . . ., the new quantum number N — —2, —1, 0, . . .. For all quantum numbers N a system of (up to 
eight) coupled differential equations has to be solved. For N = — 2 the system is reduced to one equation that 
corresponds to a state with purely heavy hole character. Non-axially symmetric systems can be treated by taking the 
coupling between the solutions of the axially symmetric problem [Eq. I|14|) ] into account>2S The form of the coupling 
depends on the symmetry along the growth direction and can be included by the substitution of a linear combination 
of the wave functions & 



J2cN^ N {r), 

N 



(15) 



with K = -2, -1, 0, 1 and N = K, K + 4, K + 8, . . . for a (001) oriented structure (C 4 symmetry); K = -2, -1, and 
TV = K, K + 3, K + 6, . . . for a (111) structure (C 3 symmetry); K = -2, -1 and N = K, K + 2, K + 4, . . . for a (110) 
structure (C2 symmetry); and K — —2 and N = K, K + 1, K + 2, . . . for other growth directions (C\ symmetry). A 

system of coupled differential equations for the envelope functions /j^' has to be solved for each value of the quantum 
number K. The Hamiltonian matrix elements (ip ni \Hij\ip nj ) are determined using Eqs. (|llf) . (|12|) . 

The Landau level spectra of a HgTe/Hg n .3Cdo.7Te (001) QW calculated with and without applying the axial ap- 
proximation are shown in Fig. [21 (a) and (b), respectively, for the structure whose corresponding subband dispersion 
is presented in Fig.^(c). As far as the Landau level fan diagrams do not show a notable difference we will use the 
axial approximation in the following. As a result of the inverted band structure the lowest Landau level of the HI 
conduction subband and the highest Landau level of the H2 valence subband cross at B 1=3 14 T. Such a behavior is 
specific for type-Ill QW's and has been examined theoretically and experimentally (see, for example, Ref . l3li) . The 
lowest HI Landau level, which corresponds to the quantum number N — —2, has purely heavy hole character, while 
the other Landau levels of the HI subband are mixed states. The H2 Landau level with N — contains both heavy 
and light states. At B > 14 T this level becomes the lowest level of the conduction band. 



C. Exchange interaction in magnetic Hgi- B Mn B Te QW's 

In the presence of a magnetic field, the sp — d exchange interaction of the s and p band electrons with the 3d 5 
electrons of Mn in Hgi„j,Mnj,Te layer influences the band structure of the QW. Such an interaction can be taken into 
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1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 1 2 3 4 5 6 7 8 9 10 1 1 12 13 14 15 



B (T) B (T) 

FIG. 2: Landau levels of the E2, HI and H2 subbands for a n-type HgTe/Hgo.3Cdo.7Te (001) QW as a function of magnetic field 
(a) with and (b) without applying the axial approximation, dw = 12.2 nm, n2DEG = 3.47 x 10 12 cm -2 , tidl = —0.34 • «.2deg, 
and T = 4.2 K. The thick line represents the chemical potential. 



account by adding an appropriate exchange term (H ex ) to the Hamiltonian [Eq. (0] in accordance with Refs-ITH and 
I2H which leads to 

H + H ex =H-J2j(r-R n )aS n , (16) 

Rn 

where u is the spin operator of the band electrons at the position r, S n is the total spin operator of the nth Mn ion 
at position R n , and J(r — R n ) is the electron-ion exchange integral. Since the electron wave function is extended, 
the spin operator S n can be replaced by the thermal average over all states of Mn moments (S z ) for a magnetic field 
in the z-direction (mean field approximation). Moreover, within the virtual crystal approximation, J(r — R n ) can be 
replaced by yJ(r — R), where y is mole fraction of Mn, and the summation is now carried out over all cation sites. 
The exchange term in Eq. (|16(l then becomes^ 

H ex = -a g {S l ,)y^2j{T-R), (17) 

R 

The average (S z ) of the z component of Mn spin in the approximation of non-interacting magnetic moments is 
determined by the empirical expression^ 

^=- s M^^k)- (18) 

where B 5 / 2 (Z) is the Brillouin function for a spin of S — 5/2; gun — 2 is the g-factor of Mn; and the effective spin 5"o 
and the effective temperature (T + Tq) account for the existence of clusters and antifcrromagnetic interaction between 
Mn ions. The values for So and T are taken from the literature^ 

The matrix elements of H ex in terms of the Bloch functions [Eqs. (QJ] have the form: 
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(19) 



with 



A=--yN /3(S z ). 



(20) 



8 



Here, No is the number of unit cells per unit volume; a and /3 are constants which describe the exchange interaction 
according to the s — d and p — d exchange integrals, respectively. Experimental values for a and (5 can be found, for 
example, in Ref. l23l 

The sp — d exchange interaction changes the spin splitting of the conduction and valence bands in a magnetic held. 
In the parabolic approximation the effective g-factor for the Tq states can be described by the following equation [cf. 
Eqs. G3J), OH and HJ]: 

* yN a(S z ) 

9 *"= 9 -^B-> (21) 

where g* is the (/-factor of the band electrons (without exchange term) . The effect of the exchange interaction on the 
Tg states can be expressed by replacing the parameter k with 

yN (3(S z ) 

Ke " = K+ ~6^B— (22) 

The influence of the sp — d exchange interaction on the band structure is obvious when we compare the Landau 
levels in Fig. 0for the non-magnetic structure with that in Fig.|3|for a magnetic Hgo.9sMno.02Te ZHgo.3Cdo.7Tc (001) 
QW. The QW width, 2DEG density and temperature are the same in both cases. The subband dispersion for the 
magnetic QW under consideration is given in Fig.Q](d). The parameters N - a = 0.4 eV, Nq ■ /3 = —0.6 cV, So = 5/2 
and To = 2.6 K are taken from the literatureAS^ Due to the exchange interaction the lowest HI Landau level with 
quantum number N = —2 (which contains pure |r 8 , —3/2) Bloch components) is bent upwards for low magnetic fields. 
In contrast to the non-magnetic case, pairs of Landau levels from the HI subband cross even at moderate magnetic 
fields. At high magnetic fields the ordering of the levels is the same as for non- magnetic structures (Fig. . Such 
behavior was also reported for n-type Hgi_j,Mnj,Te mixed crystals^ The crossing of the lowest Landau level of the 
HI subband with the N = Landau level of the H2 subband occurs at lower magnetic fields (B « 12 T) due to the 
exchange enhanced shift towards higher energy of the H2 level. 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 
B (T) 



FIG. 3: Landau levels of the E2, HI and H2 subbands for a n-type Hgo.g8Mno.02Te /Hgo.3Cdo.7Te (001) QW as a function of 
magnetic field (dw = 12.2 nm, W2deg = 3.47 x 10 12 cm -2 , tidl = —0.34 ■ 712DEG, and T=4.2 K). The thick line represents the 
chemical potential. 



In order to compare the calculations with experimental data, the density of states (DOS) at the Fermi level has to 
be calculated from the Landau level spectrum (Fig. EJ) , because experimentally the Landau level structure becomes 
visible through the magnetic field dependence of the longitudinal resistance. The Shubnikov-de Haas (SdH) oscillations 
which are observed in the experiments are directly related to changes of the DOS at the Fermi energy. Assuming a 
Gaussian broadening of the Landau levels, the DOS is given by2i 

c n V \ / 

where the summation runs over all Landau levels. T = Tq-J B / B (B n =l T) is the Landau level broadening 
parameter^ In Fig.0|the calculated DOS for the Landau level spectrum of Fig.[3]is presented together with the SdH 
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measurement of a sampled Growth and transport characterization parameters (the QW width, the 2DEG density, the 
doping profile, etc.) have been used for the band structure calculations. The broadening parameter has been chosen 
to be To=l meV. The main features such as oscillation period, beating nodes and maxima are in good agreement. 
For a more quantitative comparison the magnetic field dependence of the diffusion constant had to be taken into 
account^ 



1 2 3 4 5 6 7 




1 2 3 4 5 6 7 



B (T) 

FIG. 4: Density of states of the n-type Hgo.9sMno.02Te /Hgo.3Cdo.7Te QW at the Fermi level (thick line) compared with 
experimental SdH oscillations (thin line). 



III. CONCLUSION 

A detailed description has been presented of self-consistent band structure calculations within an eight-band k ■ p 
model in the envelope function approach with a special emphasis on type-Ill HgTe/Hgi-^Cd^Te QW structures. 
This model is an important tool for the interpretation of both opticaU&Sl and transport experiments, where strong 
spin-orbit subband splitting effects are observed. The model has been adopted to account for sp — d exchange effects 
when magnetic (Mn) ions are introduced into the structures. Self-consistently calculated band structure and DOS are 
in good agreement with experimental transport results for Hgo.9sMno.02Te /Hg .3Cdo.7Te QW's. The calculated band 
structure of Hgi-yMn^Te /Hgi-^Cd^Te QW's and the ensuing comparison with experimental data make it possible 
to understand the mutual influence of the sp — d exchange interaction and the two-dimensional confinement effects 
on the transport properties. Moreover, the effect of the QW parameters (width, doping profile, etc.) on the values of 
a, f3, So and To can now be studied by a direct comparison of experimental data and band structure calculations. 
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APPENDIX A: CORRECTIONS TO THE MATRIX ELEMENTS OF THE HAMILTONIAN FOR THE 

[kkl] GROWTH DIRECTION 

The approach of Los et alm^ can be used to carry out calculations for (kkl) oriented QW's. Since the Tq states 
as well as the coupling between Tq and T% (T^) bands are spherically symmetric, only the Bloch basis functions 
Ui(r) (i = 3 ... 8) [Eqs. (gj] have to be transformed into symmetry adapted basis functions Ui(r'). In addition, the 
coordinate system is rotated to (x' ,y' ,z') such that the z'-axis is oriented along the [kkl] growth direction. The 
corresponding terms are added to the matrix elements of Eqs. ©. The corrections which depend on k and I (h = l/k) 
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are as follows (for simplicity the new coordinates are referred as x, y, z)t 

AV = AV a + AV C , 
h 2 6 



AV r _ 



2m (h 2 + 2) 2 
h 2 6 
2m~o~(h 2 + 2) 2 



AV a = - '" 77o " ,J 2h 2 + \){ l ikl-2k zi ik z ) , (Al) 



(h 2 - 1) (,i(fc 2 _ - hV2k x {», k z }) , 



AR = Ai? a + Ai? c , 

Ai?a - f / ^ n ,J 2h 2 + l) f ,k 2 _, (A2) 



2m (h 2 + 2) 2 
ft 2 V3 
2m (/i 2 + 2) 2 

hs/2 [{h 2 + 5)k+ - (h 2 - l)fc_] {/x, fc z }) , 



■ ((2h 4 + 6/i 2 + l)nk 2 + + 2(/i 2 - l)(>fc 2 - 2fc^fc z 



A^i = AS ± = AS ± = AS a± + AS C± , 

AS a± = ^(^^^MmA}, (A3) 



AS C± = 2 ^ o (^2^2)2 ( 2 M fe2 ~ l)(Atfcf - 2fc^fc z ) + Hh 2 - l)/xfc 
- /i(/i 2 + 5)Aifc 2 z + 2V2(/i 2 -l)fc T { A i,fc z } 



The above terms are separated into axial (index a) and cubic (index c) components. It can be shown that the axial 
and non-axial approximations give the same result only for (001) and (111) oriented structures at fc|| = 0. 

APPENDIX B: INFLUENCE OF STRAIN AND PIEZOELECTRIC EFFECTS 

The effects of strain due to the lattice mismatch between HgTe and Hg 1 _ :E Cd x Te can be taken into consideration 
by applying a formalism introduced by Bir and Pikus. 35 Terms proportional to the strain tensor e are added to the 
matrix elements of the Hamiltonian [Eq. (5)]; H nn > + H^f,. The Bir-Pikus Hamiltonian H BP is derived from Eq. (5) 
by the following substitution: 

k^kj ► €ij . (Bl) 

The strain tensor components (e^ ) transform as the product kikj and are determined using the model of De Caro et 
al& The band structure parameters have to be replaced by the deformation potentials; 



2 



2m r 



H ... c 

-2a, 

-b, (B2) 

Here, C and a are the hydrostatic, and 6 and d the uniaxial deformation potentials. Due to the strain, the coupling 
matrix elements between conduction (Tq) and valence (Tg, Tj) bands have additional terms which are proportional 
to the Kane momentum matrix element These elements are actually quite small and consequently are neglected 



(2F + 


1) 


h 2 






7i 






h 2 






72 






h 2 






73 


m 
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here. The Bir-Pikus Hamiltonian for (001) oriented QW's can be written as& 



[T e 

T e 



\ 





U e + V e S e R e 
S\ U e -V e R e \ 

H = R\ U € -Ve Se \, 



R\ -S\ U e + V e 

-j-St V2V e yflSe V2R e 

v -y/2R\ J\s\ -V2V e ~^S t 



V2V ^S c 

y/%St -V2V t 

V2Rt -±S$ 

Ue 



(B3) 







/ 



where 




(B4) 



Se — d(e xz iey Z ), 

Rc ^7~b(txx ^yy) d~ idc X y- 



tr(e) = e xx + e yy + e zz is the trace of the strain tensor. 

For (kkl) oriented structures the Hamiltonian should be presented in the symmetry adapted set of basis functions 
as described in appendix A. The transformed Hamiltonian has the form of Eq. (|B3|I . with appropriate corrections to 
the matrix elements. These corrections can be derived from Eqs. (jAlfl . I|A2(1 and (|A3|) by the substitutions indicated 
in Eqs. ijBlJ) and (B^. 

If the strain tensor has non-zero off-diagonal components (shear components) , internal electric fields are generated 
in the QW due to the piezoelectric effect. We have calculated the strain- induced polarization and electric fields as 
described in Ref. [33, and have found that the influence of piezoelectric fields on the band structure of fully strained 
HgTe/Hgi-^CdxTe (112) heterostructures is negligible^ 
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